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ABSTRACT 



Context. Short supersoft X-ray source (SSS) states (durations < 100 days) of classical novae (CNe) indicate massive white dwarfs 
that are candidate progenitors of supernovae type la. 

Aims. We carry out a dedicated optical and X-ray monitoring program of CNe in the bulge of M 31. 

Methods. We discovered M31N 2007-1 la and determined its optical and X-ray light curve. We used the robotic Super-LOTIS tele- 
scope to obtain the optical data and XMM-Newton and Chandra observations to discover an X-ray counterpart to that nova. 
Results. Nova M31N 2007-1 la is a very fast CN, exhibiting a very short SSS state with an appearance time of 6-16 days after outburst 
and a turn-off time of 45-58 days after outburst. 

Conclusions. The optical and X-ray light curves of M31N 2007-1 la suggest a binary containing a white dwarf with M WD > 1.0 Mq. 
Key words. Galaxies: individual: M 31 - novae, cataclysmic variables - stars: individual: Nova M31N 2007-1 la - X-rays: galaxies 



1. Introduction 

Classical novae (CNe) are thermonuclear explosions on the sur- 
face of white dwarfs (WDs) in cataclysmic binaries that result 
from the transfer of matter from the companion star to the WD. 
The transferred hydrogen-rich matter accumulates on the surface 
of the WD until hydrogen ignition starts a thermonuclear run- 
away in the degenerate matter of the WD envelope. The resulting 
expansion of the hot envelope causes the brightness of the WD to 
rise by a typical outburst amplitude of ~ 12 magnitude s within a 
few days, and mass t o be ejected at high velocities (see lHernanzj 
120051; IWarner|[l995l, and references therein). However, a frac- 
tion of the hot envelope c an remain in steady hydrogen burning 
on the surface of the WD dStarrfield et alJll974tlSala & Hernanzj 
120051) . powering a supersoft X-ray source (SSS) that can be ob- 
served direc tly once the eje c ted envelope b ecomes sufficiently 
transparent dStarrfield 1 198ft Srautter 2002). In this paper, we 
define the term "appearance of the SSS" from an observational 
point of view as the time when the SSS becomes visible to us, 
due to the decreasing opacity of the ejected material. This is not 
the same as the actual turn-on of the H-burning on the WD sur- 
face at the beginning of the thermonuclear runaway. The turn off 
of the SSS, on the other hand, does not suffer from extinction 



* Partly based on observations with XMM-Newton, an ESA Science 
Mission with instruments and contributions directly funded by ESA 
Member States and NASA 



effects and clearly marks the actual turn off of the H-burning in 
the WD atmosphere and the disappearance of the SSS. 

The duration of the supersoft phase is related to the amount 
of H-rich matter that is not ejected and on the luminosity of 
the white dwarf. More massive WDs need to accrete less matter 
to initiate the t hermonuclear runaway, because of their higher 
surface gravity dJose & Hernanz| [l998j). In general, more mas- 
sive WDs retain less accreted mass after the explosion, although 
this also d epends on the acc retion rate, and reach higher lu- 
minosities dYaron et all 12005). Thus, the duration of the SSS 
state is inversely related to the mass of the WD. In turn, the 
transparency requirement mentioned above implies that the ap- 
pearance time is de termined by the fract i on of mass ejected 
in the outburst (see ISala & Hernanzj [20051: iTuch man & Truranl 
119981: lHachisu & Katoll2006l) . X-ray light curves therefore pro- 
vide important clues to the physics of the nova outburst, ad- 
dressing the key question of whether a WD accumulates mat- 
ter over time to become a potential progenitor for a type la 
supernova (SN-Ia). The duration of the SSS state provides the 
only direct indicator of the post-outburst envelope mass in CNe. 
For massive WDs, the expected SSS phase is very b rief (< 100 
d) (ISala & Hernanzj 120051 : ITuchman & Truranlll998|) and could 
have easily been missed in previous surveys (e .g. ISupper et al.l 
l200TtlWimams et al.ll20()4t iPietsch et al.ll2005bl) . 

There are only a few novae with very short SSS states known. 
One of them is the well-studied galactic nova RS Oph. RS Oph is 
a recurrent nova (RN), which underwent its last outburst in 2006. 
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X-ray observations documented a rapid decline in supersoft 
emission at ~ 90 days after t he optical outburst and a duration of 
the S SS phase of ~ 60 days (Osborne et al. 2006 1: lHachisu et all 
|2007|) . From theoretical models. lHachisu et al.l (120071) estimated 
a WD mass of (1.35 + 0.01)Mg. RNe are classified by the ob- 
servational fact that they have had more than one recorded out- 
burst. These objects are very good candidates for SN-Ia progen- 
itors, as RNe are believed to contain massive WDs. However, 
one of the problems with this RNe - SN-Ia connection was thei r 
rare occurrence in optical surveys dDella Valle & Livio| [i996). 
But hypothetical previous outbursts may have been missed, and 
then these objects were just classified as CNe in optical surveys. 
Interestingly, a very short SSS stage of a CN may indicate a very 
massive WD, which is capable of frequent outbursts, and thus is 
a potential RN. 

Another example for a nova with sho r t SSS state is the CN 
Nova LMC 2000, for which iGreineretall (|2003), based on non- 
detections of supersoft emissi on, suggest an SSS p hase shorter 
than seven weeks. The work of Pie tsch et al.l(l2007bl) already dis- 
cusses nova M31N 2004-1 If as an object with a very short SSS 
state of 35-55 days. This nova is an RN. 

This work presents the results of a new monitoring strat- 
egy tailored to CNe with short SSS states. We carried out com- 
bined XMM-Newton/ C/iaw ofra monitoring of the central region 
of M 3 1 (distance 780 kpc, lHollandl[T99i IStanek & Garnavichl 
1998) from November 2007 to mid-February 2008. Individual 
observations were separated by just 10 days. Section|2]provides 
detailed information on the X-ray and optical data sets. Results 
are presented in Sect.[3]and discussed in Sect.|4] 



2. Observations and data analysis 

2.1. Optical data 

The optical data used in this work were obtained with 
Super-LOTIS (Live rmore Optical Transient Imaging System, 
Williams et al. 2008), a robotic 60 cm telescope equipped with 
an E2V CCD (2kx2k) located at Steward Observatory, Kitt Peak, 
Arizona, USA. Starting in October 2007, the telescope was used 
every good night to monitor the bulge of M 3 1 . Using four Super- 
LOTIS fields (field of view: 17' x 17'), we covered an area of 
~ 34' x 34' centered on the core of M 3 1 . The pixel scale is 
0.496"/pixel and the typical limiting magnitude is 19 mag, in 
Johnson R. The data were reduced by a semi-automatic rou- 
tine and the astrometric and photometric ca libration utilizes th e 
M 3 1 part of the Local Group Survey (LGS. lMassev et al.ll2006l) . 
Typical values of Super-LOTIS astrometric and photometric lcr 
accuracies are 0'.'25 and 0.25 mag, respectively, averaged over 
the whole magnitude range. 

2.2. X-ray data 

The X-ray data were taken as part of our M 31 nova mon- 
itoring project] using the Chandra High-Resolution Camera 
Imaging Detector (HRC-I) and the XMM-Newton European 
Photon Imaging Camera (EPIC) PN detector. The individual 
observations used here are summarized in TableQ] which lists 
the telescopes and instruments used, the observation identifica- 
tions (ObsIDs), the dates, source count rates, and luminosities 
or upper limits. With the help of the mission count rate simula- 
tor PIMMS v3.9a, we computed the energy conversion factors 



1 http://www.mpe.mpg.de/~m3 lnovae/xray/index.php 



needed to transform count rates into unabsorbed fluxes for the 
individual telescopes. 

For XMM-Newton data we applied background screen- 
ing and used the XMMSAS (XMM-Newton Science Analysis 
Software) v8.0 tasks eboxdetect and emldetect to detect 
sources in the image and perform astrometry and photometry. 
For computing upper limits, we added an artificial detection at 
the position of the source to the eboxdetect list. This list was 
used as input for emldetect (with fixed positions and likelihood 
threshold of zero) to derive the observed flux, or an upper limit, 
for all objects in the list. 

We reduced the Chandra observations with the CIAO v3.4 
(Chandra Interactive Analysis of Observations) software pack- 
age. The source detection was done with the CIAO tool 
wavedetect, and an adapted version of the XMMSAS tool 
emldetect was used to estimate background-corrected and 
exposure-corrected fluxes and count rates for the detected 
sources. We applied an astrometric calibration to the object list 
for each indiv idual observation using the X-ray source catalog of 
Kaarej] d2002l) to improve the astrometric accuracy of our mea- 
surements. 



3. Results 

In our Super-LOTIS monitoring data of 2007 Nov 2.28 UT, we 
discov ered M31N 2007-1 la as a previously unknown nova can- 
didate dPietsch et al.ll2007al) . The position of M31N 2007-1 la is 
RA = 00h42m37.29s, Dec = +41°17'10'.'3 (J2000, accuracy of 
0.2"), which is 1' 19" west and 1' 2" north of the core of M 31. 
The detected R band magnitudes are 16.7 (2007 Nov 2.28), 16.8 
(2.37), 18.0 (5.27), 18.4 (5.35), 18.4 (7.35). The optical light 
curve of M31N 2007-1 lais given in Fig.Q] We estimate the pho- 
tometric accuracy of our measurements of M31N 2007-1 la to 
be ~ 0.3 mag, which is slightly less accurate than average since 
the object is located very near the center of M 31, and measure- 
ments suffer from strong background light. The outburst date of 
M31N 2007-1 la is well-constrained, since it was first detected 
on MJD = 54406.2 but not on MJD = 54405.3. Therefore, we 
estimate MJD = 54405.75 for the actual outburst, with an error 
of +0.45 days. The horizontal error bars in Fig.Q]represent this 
uncertainty for the days after the nova outburst. 

In the second HRC-I monitoring observation of our M 31 
monitoring campaign (starting 2007 Nov 17.76), CXOM31 
J004237. 3+41 1710 was detected as a relatively bright source 
(348 ± 23 counts), whereas in the first HRC-I observa- 
tion ten days earlier, no source was detected at this posi- 
tion. There is no known X-ray source within 25" of the 
position of CXOM31 J004237.3+41 1710, in previous M 31 
XMM - Newton and Cha n dra X-ray sou r ce ca ta logs dKaaret 
2002; iKong et alj 1 20021: IWilliams et all l2004t iPietsch et al. 
2005bl: IStiele et alj 120081) . The luminosity of the source de- 
clined during the three subsequent Chandra observations (Nov 
28.79, Dec 07.57, Dec 17.49). In the first (Dec 29) and all 
following XMM-Newton monitoring observations, CXOM31 
J004237.3+41 1710 is no longer detected. The X-ray light curve 
of CXOM31 J004237.3+411710 (Fig.Q) shows that this object 
is a very fast transient. Count rates and luminosities (see below) 
are given in TableQ] 

The Chandra position of CXOM31 J004237.3+41 1710 is 
RA = 00h42m37.29s, Dec = +4ri7'10'.'2 (J2000, 3<t error of 
0'.'8), which is in excellent agreement (distance of O'.'l), within 
the errors, with the optical position of M31N 2007-1 la . The X- 
ray position uncertainties are only statistical and do not include a 
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Telescope/Instrument 


ObsID 


Exp. Time n 


Start Date 


Offset* 




Count Rate 




T 

1^0.2-1 .0 






rksi 


run 






ret s _l l 




fprg s ' 1 


Chandra HRC-I 


8526 


19.9 


2007-11-07.64 


5.89 




< 8.7 xl0~ 4 




< 4.0 xlO 36 


Chandra HRC-I 


8527 


20.0 


2007-11-17.76 


16.01 


(17.4 


± 1.1) xl0~ 3 


(8.1 


±0.5) xlO 37 


Chandra HRC-I 


8528 


20.0 


2007-11-28.79 


27.04 


(4.3 


±0.5)xl0- 3 


(2.0 


±0.2)xl0 37 


Chandra HRC-I 


8529 


18.9 


2007-12-07.57 


35.82 


(4.7 


±0.5)xlO- 3 


(2.2 


±0.2) xlO 37 


Chandra HRC-I 


8530 


19.9 


2007-12-17.49 


45.29 


(3.0 


±0.4) xl0~ 3 


(1.4 


±0.2)xl0 37 


XMM-Newton EPIC PN 


0505720201 


22.2 


2007-12-29.57 


57.80 




< 8.7 xl0~ 4 




< 2.7 xlO 35 


XMM-Newton EPIC PN 


0505720301 


21.9 


2008-01-08.29 


67.53 




< 3.0 xl0~ 3 




< 5.4 xlO 35 


XMM-Newton EPIC PN 


0505720401 


18.1 


2008-01-18.63 


77.87 




< 2.2 xl0~ 3 




< 6.0 xlO 35 


XMM-Newton EPIC PN 


0505720501 


17.2 


2008-01-27.94 


87.17 




< 2.0 xl0~ 3 




< 1.3 xlO 36 


XMM-Newton EPIC PN 


0505720601 


17.3 


2008-02-07.20 


97.44 




<4.3 x!0~ 3 




<4.2 xlO 35 



Notes: Dead-time corrected 

b : Time in days after the outburst of nova M31N 2007-1 la in the optical dPietsch et al.ll2007ah on 2007 Nov 2.28 UT 

(MJD = 54406.2); uncertainties are 0.45 days, according to Sect.[3] 
c : X-ray luminosities and 3tr upper limits (unabsorbed, 0.2 - 1.0 keV) estimated using a blackbody model with kT = 50eV 

and N H = 6.7x 10 20 cirT 2 (see Sect.O 
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Fig. 1. Optical (black symbols and left y-axis) and X-ray (red 
symbols and right y-axis) light curve of M31N 2007-1 la . 
Detections are marked by crosses, the sizes of which correspond 
to the errors in magnitude/luminosity and time after outburst 
(MJD = 54405.75). Upper limits are represented by horizontal 
bars that show the uncertainty in time after outburst. Detections 
are connected with solid lines, upper limits to detections by dot- 
ted lines. 



possible syste matic offset betwee n the catalog of lKaaret) {2002) 
and the LGS dMassev et al.ll2006l) . 

In Fig.Q] we give the X-ray light curve of CXOM31 
J004237.3+41 1710 with respect to the outburst of M31N 2007- 
11a. The delay times to the nova outburst given in TableQ] for 
upper limits and detections imply that the appearance of the SSS 
happened between 6 and 16 days after the optical outburst and 
that the SSS state lasted between 29 and 52 days. 

There is no detailed spectral information about CXOM31 
J004237. 3+41 1710 , because of the very limited energy res- 
olution of the HRC detector and the source not being visi- 
ble in our XMM-Newton monitoring observations. Also, ob- 
servations with the Swift X-ray telescope from 2007 Nov 24 



to Dec 2 (ObsIDs 031028001 to 031028008) and with the 
Chandra Advanced CCD Imaging Spectrometer ACIS-I on 
2007 Nov 27 (ObsID 8187) were not deep enough to help 
with the spectral classification. (There is also no UV detec- 
tion of CXOM3 1 J004237.3+41 1710 with the Swift UV/Optical 
Telescope or the XMM-Newton Optical Monitor.) We there- 
fore used the hardness ratios from the HRC-I described in "The 
Chandra Proposers' Observatory Guide'Q for ObsID 8527 de- 
rived from the rates in the bands S, M, and H (channels 1:100, 
100:140, and 140:255, respectively). We derived ratios S/M 
= -0.10 ± 0.15, M/H = 0.09 ± 0.15, indicating a very soft 
spectrum. To compare luminosities from Chandra and XMM- 
Newton upper limits, we assumed a blackbody model with kT = 
50eV and absorption column of 6.7xl0 20 crrr 2 , the galactic 
foreground Nh towards M 31. An absorbed blackbody model 
is the simplest way of describing an SSS spectrum (see also 
iPietsch et al.l l2007bb . The temperature is an average value for 
the M 3 1 novae described in the works of IPietsch et al. (2007b, 
20051. We assume a sourc e distance of 780 kpc dHoUandll998t 



Stanek & Garnavichl fl998l) . The usage of a generic spectral 



model, Although only done to compare the different instruments, 
using a generic spectral model may of course produce signifi- 
cant luminosity uncertainties. Therefore, in Table|2]we give the 
HRC energy conversion factors (count rate / flux, computed us- 
ing PIMMS), as well as unabsorbed HRC peak luminosities and 
PN upper limits of CXOM31 J004237.3+41 1710 for different 
blackbody temperatures and an Nh fixed to the Galactic fore- 
ground value. The table shows that the source can be classified 
as a luminous SSS (Lo.2-1.0 > 10 37 erg s ) independent of the 
assumed blackbody temperature. 

4. Discussion 

From the positional and temporal agreement and from the HRC-I 
hardness ratios, we infer that CXOM31 J004237.3+41 1710 and 
M31N 2007-1 la are the same object. 

The ejected mass in the nova outburst can be estimated 
from the start date of the SSS phase and from the expansion 



see http://cxc.harvard.edu/proposer/POG/html/index.html 
chapter 7.6 
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Table 2. HRC energy conversion factors (ecf's) and unabsorbed 
HRC and PN luminosities for different blackbody temperatures. 



kT 
[eV] 



ecf™c 
[10 10 ct cm 2 erg-'] 



j HRC a 
^0.2-1.0 

[10 37 ergs" 1 ] 



I PN b 
^0.2-1.0 

[10 35 erg s- 1 ] 



40 


1.16 


10.9 ± 0.7 


3.6 


50 


1.57 


8.1 ±0.5 


2.7 


60 


2.09 


6.1 ± 0.4 


2.3 


70 


2.71 


4.7 ± 0.3 


2.1 



Notes: a : For ObsID 8527 (see Table[U 

b : For ObsID 051 1380201 (see TablefJJ 




10 15 20 

SSS appearance (days after outburst) 



Fig. 2. Relation of ejected mass (in solar units) and SSS appear- 
ance time for different expansion velocities v (see Sect. [4]). 



velocity of the ejected material. The decrease in the optical 
thickness of the expanding ejecta is responsible for the rise in 
the X-ray li ght curve of the post-outburst novae, a s shown for 
V1974 Cvg (IShore et alJll996t IKrautter et alJ [T996). Assuming 
that the material ejected by the nova explosion forms a spher- 
ical, homogeneous shell expanding at constant velocity v, the 
hydrogen mass density of the shell will evolve in time t as 
p = Mjl /(|^y 3 f 3 ) where M e ^ is the ejected hydrogen mass 
(IKrautter et al.lll996T) . The column density of hydrogen evolves 



with time as N H (cm 2 ) = l{^nuv 2 t 2 ), where u = 1.673 x 
10~ 24 g is the atomic mass unit. We assume a typical value for 
the expansion velocity of 2000 km s , since we do not have an 
optical spectrum for M31N 2007-1 la , and that the SSS turns 
on when the absorbing hydrogen column density decreases to 
~ 10 21 cm 2 . With this, the appearance time of the SSS source 
(between 6 and 16 days) constrains the ejected mass to the range 
(0.4 - 3) x 10" 7 Mq. This is about two orders of magnitude lower 
than ty pical ejected masses determined in a similar way for M3 1 
novae (Pietsch et al. I l2007bl) . However, the unknown expansion 
velocity error adds additional uncertainties to the estimate of 
M e ^. Figure|2] shows the range of M e ^ for different expansion 
velocities as a function of the SSS appearance time. 

The turn-off time of the SSS phase allows us to estimate the 
amount of hydrogen-rich material burned on the WD surface, 
M b " m = (LAt)/(X H e), where L is the bolometric luminosity, Af 
the duration of the SSS phase, Xh the hydrogen fraction of the 
burned material, and e = 5.98 x 10 18 erg g -1 . We compute the 



burned mass during the short SSS phase of M31N 2007-1 la , 
assuming a bolometric luminosity of 3 x 1O 4 L and a hydrogen 
mass fraction Xh = 0.5. The SSS turn off at 45-58 days after 
the optical outburst constrains the burned mass to the range (8 - 
10) x 1O~ 8 M0. It is noteworthy that this is within a factor 2-3 the 
same as the ejected mass. We used typical values for L and Xh, 
since we do not know the actual values for M31N 2007-1 la . 

M31N 2007-1 la has the shortest SSS phase observed so far, 
which can be caused either by a WD with Mwd > 1 ■ 1 MQfor a 
standard hydrogen fraction in the envelope or by very hydrogen- 
poor envelope in a WD with Mwd ~ 1.0 Mq. The mass of the 
stable H-burning envelope decreases with increasing WD mass, 
with massive WDs showing the shortest SSS phases. Massive 
WDs also eject less massive shells during outburst, leading to a 
faster rise in the SSS emission. 

From the optical light curve we find a decline rate of ~ 0.6 
mag/d. Extrapolating the magnitude to the day before the first 
observation (MJD = 54405.3) leads to an estimated maximum 
magnitude of m« ~ 16.1, which implies an absolute magni- 
tude at maximum of Mr ~ -8.4 ± 0.4 at the distance of M 31. 
Within the uncertainties, this value matches the maximum mag- 
nitud e versus rate of decline (MMRD) relationship obtained 
by iDella Valle & Livid (Il995l) . Assuming that the strength of 
a nova outburst (for a given temperature of the WD and ac- 
cretion rate) mostly depends on the mass of the WD (see e.g. 
lLivioll992hlbella Valle & Livioll995l;lDella Valld2002l and ref- 
erences therein), one c an estimate W D masses from optical nova 
light curves. However. iLiviol dT992l) mentions that great caution 
should be exercised when applying his maximum magnitude vs. 
WD mass relation to individual novae, since there are weaker 
dependences of the nova luminosity on the magnetic field of 
the WD and the enrichment of the envelope by heavy elements, 
which are neglected in his formula. Using our optical light curve, 
we compute thatM31N 2007-1 la likely originated from a WD 
with Mwd ~ 0.9-1.15 Mq. This result is consistent with the X- 
ray observations. 

As mentioned before, many of the novae with known short 
SSS phases turned out to be RNe. To check for the possibility 
that M31N 2007-1 la is an RN, we looked for historical novae 
with positions close to M31N 20 7-1 la. In the compiled on- 
line nova catalog of iPietsch et al. (2007bj3 there are two his- 
toric al novae with i n a r adius of ~ 20" around M31N 2007- 
11a. ISharov et alJ (|200§ fou nd a nova (M3 1N 1998-07n) at 
a distance of about 15", and iHubbld d 19291) reported a nova 
(M31N 1920-10a) at a distance of about 20" to M31N 2007- 
11a. The J2000 coordinates of M31N 1920-10a wer e taken from 
the C ombined General Catalogue of Variable Stars dSamus et al.l 
2004). From the finding chart given bv lSharovetalJ(l2000l) . nova 
M31N 1998-07n and M31N 2007-1 la do not seem to be iden- 
tical. This agrees w ith the position accuracy of 0"5 given by 
ISharov et alf (l2000). For M31N 1920- 10a no finding chart exists 
and the position given bv lHubbld (119291) has a minimum accu- 
racy of '. 1 . The actual unknown position accuracy of historical 
novae leaves the (small) possibility open that both novae in fact 
originated on the same WD. A relatively large position error for 
a historical nova could be possible, since e .g. the re-analysis of 
photographic plates by iHenze et al.l J2008) gives several exam- 
ples for historical novae with position errors of ~ 20" or more. 
Clarification of this point could only come from re-exa mination 
of the original photographic plates used in the work of Hubble 
(1929). 



http://www.mpe.mpg.de/~m3 lnovae/opt/m3 1/index.php 
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